Order-Parameter Correlation Functions in Quantum Critical Phenomena 
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We investigate the functional form of the order-parameter (two-point) correlation function in 
quantum critical phenomena. Contrary to the common lore, when there is no particle-hole symmetry 
we find that the equal-time correlation function at criticality does not display a diverging correlation 
length. We illustrate our conclusions by Monte Carlo calculations of the quantum rotor model in 
d — 2 space dimensions. 

PACS numbers: 05.70. Jk,73. 43. Nq, 05.70.Fh, 05.30.Jp 



o 

CD 

e 

I 



o 



> 

ON 

o 
p 

o 
i> 
o 



X 



To probe the dynamical response of an interacting 
quantum many-body system one typically investigates its 
response to an externally applied time-dependent force 
that perturbs the system slightly away from equilib- 
rium. The two-point time-dependent correlation func- 
tion (TPCF) is the quantity which encodes the physi- 
cal response of the system to that perturbation in lin- 
ear response theory [ij. In this way, the TPCF contains 
information about the nature of the excitations of the 
unperturbed system. If the system is close to a continu- 
ous (thermodynamic or quantum) phase transition, un- 
der certain assumptions (e.g. long-distances), the order- 
parameter TPCF may acquire a universal form indepen- 
dent of the microscopic details of the interactions. This 
is the beauty and power of the general theory of critical 
phenomena. Establishing this universal functional form, 
and determining the exponents involved, is a common 
theme in classical and quantum critical phenomena Q . 

Quantum Critical Phenomena has been introduced as 
an strairiitforward extension of Classical Critical Phe- 
nomena [35. The motivation is clear and can be traced 
back to the formal functional-integral mapping between 
a quantum problem in d space dimensions and a classical 
statistical mechanics problem in d -I- 1, where the extra 
dimension corresponds to the (imaginary) time axis r 
(see, for example, Q)- Space and time coordinates do 
not necessarily satisfy Lorentz symmetry. A measure of 
the asymmetry is quantified by the dynamic critical ex- 
ponent z which is equal to unity when Lorentz invariance 
is present. Common wisdom dictates that at the critical 
point gc continuous phase transitions display a correla- 
tion length ^ diverging as ^ \g — .gd " with critical ex- 
ponent !/ > 0, and control parameter g (temperature in 
the case of thermodynamic phase transitions). Similarly, 
the correlation time diverges as Tc ~ This fact re- 
flects the scale invariance of the system at criticality. It is 



then emphasized that consequently the equal-time (order- 
parameter) TPCF G(r, T ~ t' — 0) should decay alge- 
braically with the distance r = |r|, G(r, 0) ~ i^j.d-\-z-2+7j ^ 
thus defining a characteristic exponent rj (anomalous di- 
mension). Different relations or laws are satisfied by the 
various exponents as a result of scaling hypotheses (e.g. 
in the free energy) [sj. 

The purpose of this paper is to note that the common 
lore depicted above is not always mathematically justified 
in quantum critical phenomena and, most importantly, 
inadequate in certain cases of physical relevance. The 
root of the discrepancy lies in the lack of particle-hole 
symmetry in the case of non- vanishing chemical potential 
/i, which leads to a non-divergent length at criticality and 
short-range equal-time correlations as determined from 
G(r, 0). This in turn implies that when /.< 7^ there is no 
power law scaling ansatz in G(r, 0) at the critical point. 
On the other hand, the zero-frequency order-parameter 
correlation function, G(r, uj — 0), does display long-range 
power law behavior at criticality with a divergent corre- 
lation length ^ for arbitrary fi. This illustrates the funda- 
mental physical distinction between equal-time and zero- 
frequency correlations. We will present both a theoretical 
proof of these statements within the Ginzburg-Landau- 
Wilson (GLW) paradigm of quantum phase transitions, 
with an analytic calculation for the quantum spherical 
model, and a numerical confirmation in the quantum ro- 
tor model. Our conclusions rigorously apply as long as 
the dynamics of the order parameter of a given micro- 
scopic quantum model is faithfully described by this kind 
of GLW effective field theory. This holds regardless of 
the nature of the original degrees of freedom, i.e whether 
bosonic or fermionic. Moreover, similar derivations can 
be performed with multicomponent order parameters. 

Our starting quantum model is of the familiar phe- 
nomenological GLW form 



S[i^,r] = dr|^d'^r{-^*(r,r)p.-M)' + V2]^(r,r)+ro|V(r,T)p + |^|V'(r,r)n , (1) 
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where is the space and (EucHdean) time- 

dependent complex order parameter field, rg and uq are 
physical parameters, V = L"^ is the volume, and T is 
the temperature. The partition function is given by 
Z = J V[i;]V[ij*]exp{-S[ij,ij*]). This (Euclidean) ac- 
tion S is symmetric under the transformation ip{r,T) 
— V'(r,T), but it is not particle-hole symmetric because 
/i is, in general, not zero. This GLW functional repre- 
sents an analytic expansion (in the vicinity of the critical 
point) in '0(r, t) with interaction terms respecting the 
group of symmetries G of the disordered phase. 

Under the assumption of (bosonic) periodic boundary 
conditions along the time axis, one can Fourier transform 
the fields as 



V'(r,T) 



(2) 



with Matsubara frequencies uj = 2TrTn (n 
0,±1,±2,---), and discrete momenta ka 
2TTPa/L {pa = 0,±1,±2,---) with a = I,--- 
Then, 
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I ' dT [ d''r|V(r,T) 
Jo Jv 



(3) 



Introducing an auxiliary scalar field (j)[v, r), we decouple 
the last term by the Hubbard-Stratonovich identity 



-^/dr/d'^r|V(r,r)|'' 



const. 



with the result for the modified action 5'['0, "0*1 0] 
S = \ Y}(u:-ty)^ + k^ + r^m\^,u)\' 
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[ d'^r[-L^2(r, r) + z0(r, r)|^(r, t)^^) 
Jv 2mo 



In order to investigate the order-parameter TPCF ex- 
plicitly, we will adopt a saddle-point approximation to 
evaluate the integral over (j). The resulting (quantum 
spherical model) action is 

~S = i^[(c^-z;,)2 + fc2 + f^o]|^(k,o.)p, 



where rg = ro+i4>. For the sake of clarity, in the following 
we will concentrate in the T = and infinite volume case. 
The resulting self-consistent equation is 
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2tt J {27Ty [{uj - i^if + k^+ fo] ' 



(6) 



Equivalently, we can define a correlation length ^ 



To - At +uo 
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(7) 



where ^~ = tq + icj) — . Clearly, this integral dis- 
plays an ultraviolet divergence that needs to be regular- 
ized. To avoid this divergence, here we adopt a lattice 
regularization by replacing fc^ ^a'^i^ ^ cosfc^) and 
/ (2tt)'^ ~^ (l/^)X]k' which correctly takes into account 
the long-wavelength contributions while providing an in- 
trinsic cutoff due to the lattice constant. 

Figure [T] shows the behavior of ^ in two and three di- 
mensions, where, without loss of generality, we set uq = 1 
and use ro as the control parameter to signal the quan- 
tum phase transitions for n — and /i 7^ 0. We see 
that ^ diverges as one approaches the quantum critical 
point. Below we show that ^ is actually the correlation 
length, which gives the diverging (or singular) behavior 
at criticality, but it is not the length characterizing the 
TPCF. 

With this correlation length ^, the order-parameter 
(connected) TPCF is given by 

G(r,T) = (^(r,T)V'*(0,0))-(^(r,r)}(V'*(0,0)) 



d'^kduj 



{2TT)d+i {uj - 2/i)2 + fc2 + + ^"2 



for d > 1, where A^^ = ^2 _,_ ^-2^ p ^ (^r'^ + t2)V2/a, 
and Knip) is the modified spherical Bessel function of 
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FIG. 1: (color online) The correlation length, ^, determined 
by the self-consistent Eq. ((TJ , diverges as one approaches the 
transition point by tuning ro. Here we consider the three- 
and two-dimensional cases, (uo = 1) 
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the third kind whose asymptotic form is 



Kn{p) 



#(f) 



e >^ 
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p> 1; 
p< 1, 



(9) 



where T(n) is the Gamma function. Therefore, for fi ^ 0, 
the long-range behavior of the TPCF near the critical 
point (^~^ 0) reduces to 



G(r,T) 



5 — /i(\/''^+T^ — r) 
(^2 +^2)d/4 



(10) 



while, for /i = 0, G(r,r) - l/(r2 + r2)('i-i)/2. We note 
that, for /X ^ 0, the equal-time spatial correlation shows a 
short-range behavior characterized by a length A ^ 1///, 
not by the diverging ^, while the temporal correlation 
does display long-range behavior. 

The long-range correlation associated with the diver- 
gence of however, leads to the susceptibility thermo- 
dynamic sum rule 3| so that (^(k, w) = l/[{uj — ifi)^ + 
k"^ + 11^ + ^^2] ^ oo as ^ — > oo in the long-wavelength 
(fc — > 0) and low-frequency (w 0) limit. Therefore, the 
range of the zero-frequency correlation function is char- 
acterized by ^ in the form 



G(r,w = 0) 



exp(-r/e)/r('^-i)/2 r/^ » 1; 

,d-2 „ie ^ V-L-LJ 



1/r'' 



rli « 1, 



which, at criticality, is often compactly represented as 
G(r, a; = 0) exp(— r/f)/r''~2. Note that this form is 
the same as that of the d-dimensional classical critical- 
ity. We may view that the long-range spatial correlation 
at zero-frequency is brought by the long-range temporal 
correlation, even though the equal-time correlations are 
short ranged. 

We see that in the model a non-vanishing /z brings in 
an imaginary term linearly coupled to w in the action. 
Without this term, the quantum fluctuations along the 
temporal and spatial directions are isotropic, yielding a 
dynamical critical exponent = 1, and the mapping of 
the d-dimensional quantum phase transition to the (d -\- 
z)-dimensional classical phase transition is well justified. 
For /Lt 7^ 0, the term breaking the particle-hole symmetry 
causes an asymmetry in the temporal direction leading 
to z = 2 (i.e. there is no Lorentz invariance). This 
asymmetry property should not change at higher orders 
in uq. The formal solution, including the effect of uq, is 



G(k,cj) 



1 



{uj - ifiY -(- A;2 -I- ro - S(ro, uq, k, w) 



,(12) 



where the self-energy I](ro, uq, A^; k, w) ^ as uq — > 
and S(7'oc, uo, 0, 0) = rpc — /^^ at the critical point rg = 
roc- Therefore, we expect that the self-energy term does 
not alter the nature of the broken particle-hole symmetry 
caused by a non-zero /z, although it may introduce an 
anomalous exponent rj 0. 



Our analytic results are now confirmed by a Monte 
Carlo calculation of the quantum rotor model 



if - 2 J ^ cos( 

(r,r') 



0r'), (13) 



where U represents the rotational energy scale of a ro- 
tor at site r, J denotes the coupling strength between 
nearest-neighbor sites (r, r'), and a non-zero n{n = n/U) 
breaks the particle-hole symmetry. In the following, we 
will show that the correlations of this model near critical- 
ity are correctly captured by the effective (GLW) model 
ofEq. mi). 

We performed Monte Carlo calculations of the total 
TPCF Gt(r,r) = (V'(r, r)V'*(0, 0)) = (e^[e(r,r)-e(o,o)]^ 
a square lattice of size L x L x Lt, with periodic bound- 
ary conditions, where L is the size in a spatial direction 
and Lr in the temporal direction. For Monte Carlo cal- 
culations, we represent the quantum rotor model in the 
form of a classical {d + l)-dimensional action 



S[3] 



1 
2K 



VJ=0 

E {J. 

(r,T) 
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x{r.T) 



•Jyir.T) 



{Jr 



{r,r) 



if}, (14) 



where K ^ \J2J jV is the tuning parameter control- 
ling the quantum fluctuations, Jq,'s are integers, and 
V • J = represents the current conservation condition 
at each site. We use a recently developed worm algorithm 
to update the current configurations in Eq. ^4]) and 
measure the correlation functions. 

Figure[2]shows the TPCF as a function of a;, a distance 
along a spatial axis, in a lattice of size 40 x 40 x 400 in 
the vicinity of the critical point = 0.2978, which is 
determined by finite-size scaling of the superfluid stiff- 
ness, for \ilU = 0.20. Inspired by Eq. (fTO|) . we use the 
fitting function 



Gt(x,r = 0) = C-f a( + — -), (15) 

V a; (L — x) / 

for 1 ^ a; ^ i and obtain I « 1.4 — 1.9 for different 
if's and sizes, while G = |(V')P depends sensitively on 
K. This result strongly supports the fact that the equal- 
time TPCF has a short-range behavior, characterized by 
a finite rather than a diverging correlation length. 

To understand the nature of £, we measure it as a func- 
tion of /i at the critical point in Fig. [3l Our results clearly 
show that 1/^ is proportional to /i, and I remains finite 
except at integer fillings (/i = 0) where the particle-hole 
symmetry is restored. 

The long-range nature of one-particle correlations can 
be observed in G(r,w = 0). In Fig. [H the zero- 
frequency total correlation function, defined by Gt(a;, w = 
0) = X]r=i ^t(^, '''), displays long-range behavior near 
the critical point. We can easily see that, by subtracting 
out the constant part to obtain G(r,a; = 0), the correla- 
tion range is still beyond the size of the system. 




FIG. 2; The TPCF G{x, 0) in a lattice of size 40 x 40 x 400 for 
different T^'s near tfie critical point Kc = 0.2978 for ^/U — 
0.20. The characteristic length £, governing the range of the 
correlation function, has a finite value and varies little for 
different 7^'s and sizes. Inset: The same functions ed K = Kc 
in lattices of size 16 x 16 x 64, 20 x 20 x 100, 24 x 24 x 144, 28 x 
28 x 196, 32 X 32 X 256, and 40 x 40 x 400. 



FIG. 3: The inverse length, l/£, at the critical point 
(5^^ — 0) displays the behavior 1/1 oc n/U. Calcula- 
tions are performed in a lattice of size 40 x 40 x 400 at 
Kc = 0.3286,0.3205,0.3102,0.2978,0.2830, and 0.2655 for 
fi/U = 0.05,0.10,0.15,0.20,0.25, and 0.30, respectively. The 
dotted line is a guide to the eye. 



In summary, we have shown that, when particle-hole 
symmetry is broken, the physical quantity that displays 
a diverging length at criticality is the zero-frequency and 
not the equal-time correlation function. The latter is 
short-ranged and within our paradigm of quantum phase 
transitions it is not Lorentz invariant. We have not only 
proved this statement analytically but also illustrated its 
correctness by studying the planar quantum rotor model. 
This model shares the same critical properties as the 
Bose-Hubbard model [1] assuming that only phase fluc- 
tuations are responsible for its critical behavior. This 
observation is quite relevant given the current interest 
in the experimental simulation of quantum phase transi- 
tions in optical lattices of ultracold atomic gases [H. 

M.-C.C. thanks the hospitality of the Department of 
Physics at Indiana University where parts of this work 
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Fund Grant No. KRF-2005-041-I01625. 
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FIG. 4: Zero-frequency correlation function near the critical 
point. Long-range correlations are clearly observed. 
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